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1 Introduction and main results 



The purpose of the present paper is to put into evidence a new formula in Hanultonian dynamics, both simple 
and general, relating the time evolution of localisation observables to the variation of energy along classical 
orbits. 

Our result is the following. Let M be a (finite or infinite-dimensional) symplectic manifold with symplectic 
2-form (jj and Poisson bracket { • , • }. Let H e (M) be an Hamiltonian on M with complete flow {(pt}tm- 
Let $ = ($1 , . . . , $d) e C°° (M; M**) be a family of observables satisfying the condition 

{{^j,H},H]=Q (1.1) 

foreachj e {1, . . . Then we have (see Theorem 3.3, Corollary 3.4 and Lemma 3.6 for a precise statement): 

Theorem 1.1. Let H and $ be as above. Let f ^ Cbe such that f = 1 ona neighbourhood of 0, f = 

at infinity, and f{x) = ,f{—x) for each x G R'^. Then there exist a closed subset Crit(_ff, $) C M and an 
observable Tf G C~ (M \ Crit(if, $)) satisfying {Tf,H} = lonM\ Cnt{H, $) such that 

POO 

lim i / dt [(/($/r) o cp_t) (m) - (/($/r) o ^t) (m)] = T/(m) (1.2) 

for each m e M \ Crit(iJ, $). 

The observable T/ admits a very simple expression given in terms of the Poisson brackets djH := {^j,H} 
and the vector V-ff := {diH, . . . , ddH), namely, 

Tj :=-$.(Vi?/)(Vi/), (1.3) 

where Vi?/ •.R'^^C'^ is some explicit function (see Section 2). 

In order to give an interpretation of Formula (1.2), consider for a moment the situation where M := 
T*M." ~ R^" is the standard symplectic manifold with canonical coordinates {q,p) and 2-form w := dg-' A 

dpj . Furthermore, let H{q,p) := h{p)bca purely kinetic energy Hamiltonian and let ^{q,p) := q be the stan- 
dard family of position observables. In such a case, the condition (1 . 1) is readily verified, the vector VH reduces 
to the usual velocity observable V/i associated to H, and the l.h.s. of Formula (1.2) has the following meaning: 
For r > and m £ M \ Crit(i?, fixed, it is equal to the difference of times spent by the classical orbit 
{97((m)}(gR in the past(firstterm)and in the future(secondterm)within the region Sr := supp[/($/r)] C M 
defined by the localisation observable /(<l>/r). Moreover, if we interpret the map ^ := {Tj, •} as a derivation 
on C°° (M \ Crit(-ff, $)), then Tf on the r.h.s. of (1.2) can be seen as an observable "derivative with respect to 
the energy H" on M \ Crit(fl', since -^(H) = {Tj, H} = lonM\ Cnt{H, Therefore, Formula (1.2) 
provides a new relation between sojourn times and variation of energy along classical orbits. Evidently, this 
interpretation remains valid in the general case provided that we consider the observables as the components 
of an abstract position observable $ (see Remark 3.7). 

Our interest in this issue has been aroused by a recent paper [31], where the authors establish a similar 
formula in the framework of quantum (Hilbertian) theory. In that reference, -ff is a selfadjoint operator in a 
Hilbert space H, ^ = ($i, . . . , ^d) is a family of mutually commuting selfadjoint operators in H, (1.1) is a 
suitable version of the commutation relation , iJ ] , iJ ] =0, and Tf is a time operator for H (i.e. a symmetric 
operator satisfying the canonical commutation relation [Tf, H] = i). So, apart from its intrinsic interest, the 
present paper provides also a new example of result vaUd both in quantum and classical mechanics. Points of 
the symplectic manifold correspond to vectors of the Hilbert space, complete Hamiltonian flows correspond to 
one-parameter unitary groups, Poisson brackets correspond to commutators of operators, etc. (see [1, Sec. 5.4] 
and [24] for the interconnections between classical and quantum mechanics). Accordingly, we try put into 
light throughout all of the paper the relation between both theories. For instance, we link in Remark 3.5 the 
confinement (resp. the non-periodicity) of the classical orbits {(pt(m)}t£K, rn e M, to the affiliation of the 
corresponding quantum orbits {e**^ ^}teK, ^/^ € 7^, to the singular (resp. absolutely continuous) subspace of 
%. Moreover, we show in Section 4.5.2 that the Hilbertian space theory of [31] can be recast into the present 
framework of symplectic geometry by using expectation values. 
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We also mention that Formula (1.2), with r.h.s. defined by (1.3), provides a crucial preliminary step for the 
proof of the existence of classical time delay for abstract scattering pairs {H, H + V} (see [10], [14, Sec. 4.1], 
and [35, Sec. 3.4] for an account on classical time delay). If V is an appropriate perturbation of H and S is the 
associated scattering map, then the classical time delay t(to) for m e M defined in terms of the localisation 
operators /($/r) should be reexpressed as follows: it is equal to the l.h.s. of (1 .2) minus the same quantity with 
m replaced by S{m). Therefore, if m and 5(m) are elements of M \ Crit(iJ, then the classical time delay 
for the scattering pair {H, H + V} should satisfy the equation 

T(m) = {Tf -Tfo S){m). 

Now, the property {Tf, H}{m) = 1 imphes that Tf{m) = {Tf o ipf)[m) — t for each t e M. Since S commutes 
with ift, this would imply that 

T(m) = [{Tf -TfoS)o ^t] (m) 

for all t gM., meaning that the classical time delay is a first integral of the free motion. This property corresponds 
in the quantum case to the fact that the time delay operator is decomposable in the spectral representation of the 
free Hamiltonian (see [32, Rem. 4.4]). 

Let us now describe more precisely the content of this paper In Section 2 we recall some definitions 
in relation with the function / that appear in Theorem 1.1. The function Rf is introduced and some of its 
properties are presented. Then we prove various versions of Formula (1.2) in the particular case where the 
functions $ o (pj-^ : M ^ M** are fixed vectors x±ty,x,y gR''' (see Proposition 2.3, Lemma 2.4 and Corollary 
2.6). 

In Section 3.1, we introduce the Hamiltonian system (M, w, H) and the abstract position observable 
Then we define the (closed) set of critical points Crit(iJ, associated to H and $ as (see [31, Def. 2.5] for the 
quantum analogue): 

Cr\t{H, $) := {to e M I (Vi/)(m) ^ O}. 

When H{q,p) = h{p) and ^{q,p) := g on M = R^", Qr\t{H,^) coincides with the usual set Crit(if) of 
critical points of the Hamiltonian vector field Xh, i.e. 

Crit(ff) = {to e M I Xaim) = O} = {{q,p) G M^n | (v/i)(p) = O} = Crit(ff,$). 

But, in general, we simply have the inclusion Crit(iJ) C Cr\t{H, 

In Section 3.2, we prove the main results of this paper. Namely, we show Formula (1.2) when the localisa- 
tion function / is regular (Theorem 3.3) or equal to a characteristic function (Corollary 3.4). We also establish in 
Theorem 3.8 a discrete-time version of Formula (1.2). The interpretation of these results is discussed in Remarks 
3.5 and 3.7. 

In Section 4, we show that our results apply to many Hamiltonian systems (M, to, H) appearing in lit- 
erature. In the case of finite-dimensional manifolds, we treat, among other examples. Stark Hamiltonians, ho- 
mogeneous Hamiltonians, purely kinetic Hamiltonians, the repulsive harmonic potential, the simple pendulum, 
central force systems, the Poincare ball model and covering manifolds. In the case of infinite-dimensional man- 
ifolds, we discuss separately classical and quantum Hamiltonians systems. In the classical case, we treat the 
wave equation, the nonlinear Schrodinger equation and the Korteweg-de Vries equation. In the quantum case, 
we explain how to recast into our framework the (Hilbertian) examples of [31, Sec. 7], and we also treat an 
example of Laplacian on trees and complete Fock spaces. In all these cases, we are able to exhibit a family 
of position observables $ satisfying our assumptions. The diversity of the examples covered by our theory, 
together with the existence of a quantum analogue [31], make us strongly believe that Formula (1.2) is of nat- 
ural character. Moreover it also suggests that the existence of time delay is a very conmion feature of classical 
scattering theory. 

2 Integral formula 

In this section, we prove an integral formula and a summation formula for functions on W^. For this, we start 
by recalling some properties of a class of averaged localisation functions which appears naturally when deahng 



4 



with quantum scattering theory. These functions, which are denoted Rf, are constructed in terms of functions 
/ S L°°(M'*) of localisation around the origin eRVrhey were already used, in one form or another, in 
[17, 31, 32, 37, 38]. We use the notation {x) := ^/l + \x^ for any a; e W^. 

Assumption 2.1. The function f G L°°(IR'^) satisfies the following conditions: 

( i) There exists p> such that \f{x)\ < Const, {x) ~i' for almost every x G M**. 

( ii) f = lona neighbourhood of 0. 

It is clear that liuir^^ f{x/r) = 1 for each a; e R"^ if / satisfies Assumption 2.1. Furthermore, one has 
for each a; e M** \ {0} 



r ^ [f{„x) - X[o,i](m)] < /' ^ \f{l^x) - 1| + Const, d/x 
Jo M Jo Ji 



IJL-^^+P^ < 00, 



where x[o,i] denotes the characteristic function for the interval [0,1]. Therefore the function i?/ : ]R''\{0} C 
given by 



Rf{x) := r°°^[/(^x)-X[o,i](/x)] 
Jo M 



is well-defined. 

In the next lemma we recall some differentiabiUty and homogeneity properties of Rf. We also give the 
explicit form of Rf when / is a radial function. The reader is referred to [38, Sec. 2] for proofs and details. The 
symbol y^K^) stands for the Schwartz space on W^. 

Lemma 2.2. Let f satisfy Assumption 2.1. 

(a) Assume that §ir{x) exists for all j € {1, . . . , rf} and x € M^, and suppose that there exists some p > 
such that I ^ (a;) | < Const, (a;) for each x e M"^. Then Rf is dijferentiable on \ {0}, and its 
gradient is given by 



/•oo 

{VRf){x)= / dp{Vf){,ix). 
Jo 

In particular, iff e yiW^) then Rf belongs to C°°(M'* \ {0}). 



(b) Assume that Rf belongs to C^{W^ \ {0})/or some m > 1. Then one has for each x \ {0} and 
t > Othe homogeneity properties 

x-{WRf)ix)=-l, (2.1) 
t\"\{d"Rf){tx) = {d^Rf){x), 

where a e N*^ is a multi-index with 1 < |q;| < m. 

(c) Assume that f is radial, i.e. there exists fo G L°°(R) such that f{x) = /o(|a;|) for almost every x G M."^. 
Then Rf belongs to C°°{R'^ \ {0}), and {WRf){x) = -x-'^x. 

In the sequel, we say that a function / : M'' — )• C is even if f{x) = f{—x) for almost every x G R'^. 

Proposition 2.3. Let f -.W^ ^ Cbe an even function as in Lemma 2.2. (a). Then we have for each x gR'^ and 
each yGR'^\ {0} 



lim 
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f^t[f(^)-f{^)]=-x.iVRfM. (2.2) 



In particular, iff is radial, the l.h.s. is independent of f and equal to {x ■ y)/y^ 
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Proof. The change of variables :— t/r, u := 1/r, and the fact thai / is even, gives 

pOO 

liml/ dt[/(^)-/(^)] 

/>oo 

= lim i / dii{l[f{vx - iiy) - f{-^iy)\ - i [f{vx + ^ly) - f{^iy)\ }. (2.3) 
By using the mean value theorem and the assumptions of Lemma 2. 2. (a), one obtains that 



l\f{vx ± ny) - /(±/xy)| < Const, sup {(,ux ± fxy) ^^^''^ 

«e[o,i] 

for some p > 0. Therefore, if /x is big enough, the integrant in (2.3) is bounded by 

Const.(/x|y|- |a;|)"^^+''^ 

for all 1/ e (0, 1). This implies that the integrant in (2.3) is bounded uniformly in z/ e (0, 1) by a function 
belonging to ([0, oo), d/i). So, we can apply Lebesgue's dominated convergence theorem to interchange the 
limit on v with the integration over ji in (2.3). This, together with the fact that (V/)(— a;) = — (V/)(a:), leads 
to the desired result: 

/■oo /"OO 

lim i / At [/(^) - /(^)] = i / d/x • (V/)(-/zy) - x ■ (V/)(My)] 

/■oo 

= -/ d^ix■{Wf){^iy) 

Jo 

= -X ■ {WRf){y). □ 

The result of Proposition 2.3 can be extended to less regular functions / : M'' — )• C. The interested reader 
can check that the result holds for functions / admitting a weak derivative /' such that, for every real line 
L c R'^, /' is of class on L (see [41, Thm. 2.1.6]). We only present here the case (of particular interest for 
the theory of classical time delay) where / is the characteristic function xi for the unit ball Bi := {a; e M*^ | 
\x\ < !}• 

Lemma 2.4. One has for each a; e M*^ and each y gR'^\ {0} 



lim i / dt 



(x-ty\ /x + ty\ 



X ■ y 

y2 



r 

Proof. Direct calculations and the change of variables /x := t/r, v := 1/r, give 

/•oo /■oo 

— ^ X[-a(!^,a;,2/),y-2-a(i/,x,j/)] \^(M ± -p-j 

with a(i/, a;, y) := ^ {x^y^ — {x ■ y)^). Now, a(i/, x, y) > 0, and — a(i/, x,y) > if i/ > is small enough. 
So, the last expression is equal to 

1 / 7 r- X ■ W 

= -Vy - aW,x,y) T 
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if y is small enough. This implies that 



/•oo 
J 



1/1 / — 5 7 T X ■ y 1 / ^ ;^ r- X ■ y\ 

y2 



□ 



For the next corollary, we need the following version of the Poisson summation formula (see [15, Thm. 5] 
or [39, Thm. 45]). 

Lemma 2.5. Let g : (0, cxd) —> C be a continuous function of bounded variation in (0, cxd). Suppose that 
limt_>.oo g{t) = and that the improper Riemann integral dtg{t) exists. Then we have the identity 



/•oo /"OO 

5(0) + ^ 5(n) = / dtg{t) + 2^ dt cos{2Tmt)g{t). 

n>l Jo n>l-^0 



Corollary 2.6. Let / : R'' — > C /je an even function such that 

( i) f = 1 on a neighbourhood of 0. 

(ii) For each a G N"^ with \a\ < 2, the derivative f exists and satisfies |(3"/)(x)| < Const. (a;)~'-^+'') 
for some p> and all x G M*^. 

Then we have for each x gR'^ and each y e R*^ \ {0} 

n>l 

In particular, iff is radial, the l.h.s. is independent of f and equal to {x ■ y)/y^. 
Proof. For r > given, the function 

5.:(0,oo)^C, t^grit):=f{^)-f{^), 
satisfies all the hypotheses of Lemma 2.5. Thus 



lim i J2 [/(^) - /(^)] = i / dtgrit) + lim di cos(27rni)5. 

n>l ^ n>l ^ 



{t). 



The first term is equal to —x ■ {'VRf ){y) due to Proposition 2.3. For the second term, the change of variables 
fj, := t/r, v := 1/r, and two integrations by parts give 

/•oo 

lim / dt cos{2iTnt)gr{t) 

= lim V] / — cos{2Tm^i/v) [f{vx - fxy) - f{vx + iiy)] 

= E ii- E (^(.. - ^,) + ^(.. + ,y)^ 

J n>l \ J J / 

Ev-^ 2j/ df , 
yj lim > -- — -77— — (ux) 

j n>l ^ ' ■> 

V- V- I""" , vcos{2T:npilv) ( (92/ , 9^ . x\ 

- E l^^o E X {2.nY i ^'^^ - + + V • 
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Since ^^^-^ < oo, one sees directly that the first term is equal to zero. Using the fact that | qx^Lx (•^) I — 

Const, {x) for some p > and all x e R'', one also obtains that the second term is equal to zero. 

Therefore, 

n>l 

and the claim is proved. □ 



3 Hamiltonian dynamics 

In the sequel, we require the presence of a symplectic structure in order to speak of Hamiltonian dynamics. 
However our results still hold if one is only given a Poisson structure. A lack of examples and some complica- 
tions in infinite dimension regarding the identification of vector fields with derivations have led us to restrict the 
discussion to the symplectic case for the sake of clarity. 

3.1 Critical points 

Let M be a symplectic manifold, i.e. a smooth manifold endowed with a closed two-form w such that the 
morphism TM 9 X i-> ^H^) •= '^^^^ ^ isomorphism. In infinite dimension, such a manifold is said to 
be a strong symplectic manifold (in opposition to a weak symplectic manifold, when the above map is only 
injective; see [2, Sec. 8.1]). When the dimension is finite, the dimension must be even, say equal to 2n, and 
the 2n-form a;" must be a volume form. The Poisson bracket is defined as follows: for each / S C°° (M) we 
define the vector field X/ := (w^)-i(d/), i.e. d/( • ) = • ), and set {j,g\ := uj{Xf,Xg) for each 

/,5GC°°(M). 

In the sequel, the function H £ C°°{M) is an Hamiltonian with complete vector field Xh. So, the flow 
{(ft} associated to H is defined for alH e M, it preserves the Poisson bracket: 

{f °^t,g°¥>t} = {f,g}o'Pt, teR, 

and satisfies the usual evolution equation: 

^Jo(pt = {f,H}o<pt, tern. (3.1) 

In particular, the Hamiltonian H is preserved along its flow, i.e. H o (pt = if for all f G R. We also consider an 
abstract family $ = ($i, . . . , G C°°(M; W^) of observables\ and define the associated functions 

a^i? := G C°°(M) and VH := {d^H, . . . ,ddH) e C°°{M-M.'^). 

Then, one can introduce a natural set of critical points: 

Deiuiition 3.1 (Critical points). The set 

CuX{H, $) := {VH) 

is called the set of critical points associated to H and 

The set Crit(-ff, $) is closed in M since VH is continuous. Furthermore, since H} = d^j{XH), the 

set 

Cnt{H) := {m G M I X//(m) = 0} = {m G M | dHm = 0} 

of usual critical points of H satisfies the inclusion Crit(ff ) C Cnt{H, 
Our main assumption is the foUowing: 

need be, the results of this article can be extended to the case where H and $j are functions of class with {^j,H} also C^. 



-'({0}) C M 
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Assumption 3.2. One has {{^j,H}, _ff } = for each j <E {1, . . . , d}. 

Assumption 3.2 imposes that all the brackets {<^j,H} are first integrals of the motion given by H. When 
M is a symplectic manifold of dimension 2n, these first integrals are functions offce {l,2,...,2n — 1} 
independent first integrals Ji = H, J2, . . . , Jk (Ji, ■ ■ ■ ^ Jk are independent in the sense that their differential are 
linearly independent at each point of M)^. So, one should have {^j,H} = gj{Ji, . . . , Jk) for some functions 
Qj e C°°(M"; R). Using the properties of { • , H} as a derivation, one iirfers that {gj{Ji, . . . , Jk)~^^j, H} = 1 
outside gj{Ji, . . . , Jfe)~^({0}). Thus, if k first integrals as Ji, . . . , Jfc are known, finding functions $j satisfying 
Assumption 3.2 is to some extent equivalent to finding functions $0 solving {c&o, H} = 1 (the equivalence is 
not complete because these functions $0 are in general not since { • , H} is necessarily on Cnt{H)). 

For further use, we define the C°° -function Tf:M\ Cr\t{H, <&)-)• M by 

Tf := • {VRf){VH). 
When / is radial, Tf is independent of / and equal to 

VH 



T:=<P- 



due to Lemma 2.2.(c). In fact, the closed subset ({0}) of Af \ Crit(i7, $) admits an interesting interpretation: 
If we consider the observables $j as the components of an abstract position observable then WH can be seen 
as the velocity vector for the Hamiltonian H, and the condition 

T{m)^Q ^ $(m) • (Vi?)(TO) = (3.2) 

means that the position and velocity vectors are orthogonal at m G T~ ^ ({0}) . Alternatively, one has T{m) = 
if and only if the vector fields X|$|2 and Xh are w-orthogonal at m, that is, Wm(X|$|2(m),XH(m)) = 0. 
The simplest example illustrating the condition (3.2) is when ^{q,p) := q and H{q,p) := are the usual 
position and kinetic energy on {M, w) := (R^", Y^^=i '^1'' ^ '^Pj) ■ ^uch a case, (3.2) reduces to q ■ p = 0. 

3.2 Sojourn times of classical orbits 

Next Theorem is our main result. We refer to Remark 3.7 below for its interpretation. 

Theorem 3.3. Let H and $ satisfy Assumption 3.2. Let f -.W^ ^ Cbe an even function as in Lemma 2. 2. (a). 
Then we have for each point m G M \ Cnt{H, $) 

/oo 
[(/($/r)oy,_,)(m)- (/($/r)o(^,)(m)] =Tf{m). (3.3) 



r—^OG 







In particular, iff is radial, the l.h.s. is independent off and equal to $(m) • (yH){m)^ • 
Proof. Equation (3.1) implies that 

for each t e M. Similarly, using Assumption 3.2, one gets that 

^ H}oipt = {{$„ H}, H]o^, = 0. 
So, $j varies linearly in t along the flow of Xh, and one gets for any m G M 

o ipt){m) = {^j o ^o)(m) + 1 (-| o Vt){m)\^^J - $j(m) + t{djH){m). 



^In the setup of Liouville's theorem [5, Sec. 49], we have k = n and the first integrals are mutually in involution. Furthermore, on the 
connected components of submanifolds given by fixing the values of these n integrals in involution, the flow is conjugate to a translation 
flow on cyhnders W^'^ x (see [1, Thm. 5.2.24]). 
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This, together with Formula (2.2), gives 

^lim i / dt [(/($/r) o ^_,) (m) - (/($/r) o (m)] 

J 

= lim i n dt ^^ i'(m)+i(vg)(m) y 

= r/(m). □ 

Due to Lemma 2.4, the proof of Theorem 3.3 also works in the case / = xi- So, we have the following 
corollary. 

Corollary 3.4. Let H and $ satisfy Assumption 3.2. Then we have for each point m G M \ Cnt{H, $) 

poo 

^lim i / dt [{xii^/r) o ^_t)(m) - (xi($/r) o ^,)(m)] = $(m) • ffff/g . (3.4) 

J 

We know from the proof of Theorem 3.3 that 

{^j o ipt){m) = $j(m) + t{djH){m) for all i e M and all m e M. (3.5) 

Therefore, the l.h.s. of (3.3) and (3.4) are zero if m e Crit(i?, $). 

For the next remark, we recall that any self adjoint operator A in a Hilbert space %, with spectral measure 
■ ), is reduced by an orthogonal decomposition [40, Sec. 7.4J 

H = Hac(^) © n^iA) © -HA A) = H^c{A) ® Hs{A), 

where 'Hac(^), 'Hp(A), 'Hsc(^) and 'Hs{A) are respectively the absolutely continuous, the pure point, the sin- 
gular continuous and the singular subspaces of A. Furthermore, a vector ip G%\s said to have spectral support 
with respect to A in a set J C R if = E'"* {J)^- 

Remark 3.5. Ifm G Crit(_ff, $), then one must have ipi{'m) E Cnt{H, ^) for all f G R, since (3.5) implies 
{djH){ipt(rn)) = {djH)(rn) for all t G M. Conversely, ifm G M\Crit(if, then one must have ipt{m) ^ m 
for all t ^ 0, since $ cannot take two different values at a same point So, under Assumption 3.2, each orbit 
{(pt{iTi)}teR either stays in Cnt{H, ^)ifmG Cnt{H, or stays outside Cnt{H, $) and is not periodic if 
m i Crit(ff,$). 

In the corresponding Hilbertian framework [31 ], the Hamiltonian H and the functions $j are selfadjoint 

operators in a Hilbert space %, and the critical set n associated to H and $ is a closed subset of the spectrum of 
H. Outside k, the spectrum of H is purely absolutely continuous [31, Thm. 3.6. (a)]. Therefore, vectors € H 
having spectral support with respect to H in k belong to the singular subspace TisiH) of H, and thus lead 
to orbits {c''^ V-'}ieR confined in 'Hg(ff ) (for instance, c**^ ip stays in a one-dimensional subspace ofH ifip 
is an eigenvector of H). Conversely, vectors tp £ T-L having spectral support outside k belong to the absolute 
continuous subspace y,a,c{H) ofH, and thus lead to orbits {e**^ V'lteR contained in y,ac{H) (see [3, Prop. 5.7] 
for the escape properties of such orbits). These properties are the quantum counterparts of the confinement to 
Crit(iJ, $) (when m G Crit(ff, $)) and the non-periodicity outside Crit(iJ, $) (when m ^ Crit(-ff, $)) of the 
classical orbits {ipt{m)}tes.- 

Lemma 3.6. IfH,^ and f satisfy the assumptions of Theorem 3.3, then we have 

{Tf,H}o^, = ^^{TfOipt) = l (3.6) 

on M\ Cnt(H, $). In particular, one has Tf o ipt ^ Tf -\- 1 on M \ Crit(iJ, $). 

If we interpret the map ^ := {Tf, •} as a derivation on C°°(M \ Crit(-ff, i>)), this impUes that T/ can 
be seen as an observable "derivative with respect to the energy H" on M\ Cnt{H, since 

^(H) = {Tf,H} = l 

on each orbit {</5f (TO)}fgR, with m £ M \ Cnt{H, 
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Proof of Lemma 3.6. The first equality in (3.6) follows from (3.1). For the second one, we use successively the 
fact that ipt leaves invariant H and the Poisson bracket, Assumption 3.2, and Equation (2.1). Doing so, we get 
on M \ Crit(il, the following equalities 

^ (Tj o = - A ($ o . (Vi?/) ({$ o <p„ F}) = - A ($ + t{VH)) ■ [VRf) ({$ + t{VH), H}) 

= -^(* + t(Vif))-(Vi?/)(Vff) 
= -(Vi?) • {VRf){VH) 

= 1. □ 

Remark 3.7. Theorem 3.3 relates the sojourn times of classical orbits within expanding regions of M to the 
observable Tf. If we consider the observables $j as the components of an abstract position observable $, then 
the l.h.s. of Formula (3.3) has the following meaning: For r > and m G M \ Crit(iJ, $) fixed, it can be 
interpreted as the difference of times spent by the classical orbit {<y9t(m)}tgii in the past (first term) and in 
the future (second term) within the region '■= supp[/(<&/r)] C M defined by the localisation observable 
f{^/r). Thus, Formula (3.3) shows that this difference of times tends a.v r — > oo to the value of the observable 
Tf at m. Since Tf can be interpreted as an observable derivative with respect to the energy H, Formula (3.3) 
provides a new relation between sojourn times and variation of energy along classical orbits. 

As a final result, we give a discrete-time counterpart of Theorem 3.3, which could be of some interest in 
the context of approximation of symplectomorphisms by time-1 maps of Hamiltonians flows (see e.g. [7], [18, 
Appendix BJ, [23] and references therein). 

Theorem 3.8. Let H and $ satisfy Assumption 3.2. Let f ■.R'^ ^ Cbean even function such that 

(i) f = 1 ona neighbourhood of 0. 

(ii) For each a € N"^ with \a\ < 2, the derivative d'^f exists and satisfies \{d°'f){x)\ < Const, (x) ^^'^''^ 
for some p> and all x G M"^. 

Then we have for each point m G M \ Cnt{H, $) 

,1™ 5 E [{fi^M ° ^-n){m) - (/($/r) o <p„)(m)] = Tf{m). 

n>l 

In particular, iff is radial, the l.h.s. is independent of f and equal to $(m) • lyH){m)^ ■ 
Proof Let m e M \ Crit(il, Then we have by Equation (3.5) 

\ E ° 'P-n) (m) - (/($/r) o (m)] 

n>l 

' ^(■m)+n(V H)(my 



n>l 

and the claim follows by Formula (2.4). □ 



4 Examples 

In this section we show that Assumption 3.2 is satisfied in various situations. In these situations all the results 
of Section 3 such as Theorem 3.3 or Formula (3.6) hold. Some of the examples presented here are the classical 
counterparts of examples discussed in [31, Sec. 7] in the context of Hilbertian theory. 

The configuration space of the system under consideration will sometimes be M", and the corresponding 
symplectic manifold M = T*M" ~ M^". In that case, we use the notation {q,p), with q = {q^,. .. 
and p = {pi, . . . ,pn), for the canonical coordinates on M, and set a; := Yl^=i ^9"' A dpj for the canonical 
symplectic form. We always assume that / = Xi or that / satisfies the hypotheses of Theorem 3.3. 
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4.1 VH = g{H) 

Suppose that there exists a function g = {gi,...,gd) G C°°(R;]R'') such that VH = g{H). Then H and $ 
satisfy Assumption 3.2 since {gj{H),H} = Oforeach j. Furthemore, onehas Crit(if,$) = {g o H)~^{{0}), 
and T/ = • (V-R/) (^(-ff )) on M \ Cnt{H, $). We distinguish vanous cases: 

(A) Suppose that g is constant, i.e. g = v e R"^ \ {0}. Then Cnt{H) = Crit(iJ, $) = 0, and we have the 
equality Tf • (V-R/)(w) on the whole of M. 

Typical examples of functions H and <E> fitting into this construction are Friedrichs-type Hamiltonians and 
position functions. For illustration, we mention the case (with d = n) of H{q,p) := v ■ p + V{q) and 
<P{q,p) —qonM := M^n^ ^j^j^ ^ g M" \ {0} and F e C°°(M"; R). In such a case, onehas Vil = v 
and 

^t{q,p) = {vt + q,p- jlds{W){vs + q)). 

Stark-type Hiuniltonians and momentum functions also fit into the construction, i.e. H{q,p) := h{p)+v-q 
and ^{q,p) := p on M := M^n^ t; e M" \ {0} and /i € C°°(M";M). In such a case, one has 
Vi? = —I) and 

Ml^P) = {1 + Io^s{Vh){p-vs),p-vt). 

Note that these two examples are interesting since the Hamiltonians H contain not only a kinetic part, but 
also a potential perturbation. 

(B) Suppose that $ has only one component (d = 1), and assume that (/(A) = A for all A e M (in the Hilbertian 
framework, one says in such a case that H is ^-homogeneous [9]). Then Cnt{H, <i>) = H~^{{0}) and 
we have the equality Tf = -$(V-R/) (.ff ) on M \ H~'^{{0}). We present a general class of pairs {H, $) 
satisfying these assumptions: 

The Hamiltonian flow of the function D{q, p) := q-pon M^" is given by (f^ {q, p) = (e* q, e~* p). So, D 
is the generator of a dilations group on M^" (in the Hilbertian framework, the corresponding operator is the 
usual generator of dilations on [? (M"), see e.g. [4, Sec. 1 .2]). Therefore, the relation {£>, H} oc H holds 
for a large class of homogeneous functions H on M-^", due to Euler's homogeneous function theorem. Let 
us consider an explicit situation. Take a > and let M be some open subset of (M" \ {0}) x R". Define 
on M the function $ := and the Hamiltonian H{q,p) := h{p) + V{q), where h e C°°(R";R) is 
positive homogeneous of degree a and V G C°°(Ili" \ {0}; M) is positive homogeneous of degree —a. 
Then one has VH = {$, H} = H on M, and 

Crit(iJ) = {{q,p) e M I {Vh)ip) = (Vt/)(<z) = 0} 

C {iq,p) e M \p-{Vh)ip) = q ■ iVV)iq) = O} 
= {{q,p)eM\H{q,p)=0} 
= Crit(iJ,$). 

Furthermore, if the functions h and V and the subset M are well chosen, the Hamiltonian vector field Xh 
of H is complete. For instance, 

(i) IfV = 0, then one can take M = R^", and one has ipt{q,p) = {q + t{\'h){p),p) and 

Crit(ff) = e M I (V/i)(p) = 0} C e M I p • {Vh){p) = O} = Crit(F, $) 

(when h{p) — is the classical kinetic energy, one has Crit(if) = Crit(iJ, $) = R" x {0}). 

(ii) Let K>0. Then the Hamiltonian given by ^"(9,35) := -|- K\q\-'^) on M := R" \ {0} x R" 
has a complete Hamiltonian vector field Xh. To see it, we use the push-forward of Xh by the 
diffeomorphism (, : R" \ {0} x M" ^ R" \ {0} x R", {q,p) = (r,p), namely. 



[i4XH)]{r,p) = Y^(^{\r\'pj -2{p.ry) 



d 



(r,p) dpj 



(r,p) 
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Then, we obtain that is complete by using the criterion [1, Prop. 2.1.20] with the proper 

function gi : M" \ {0} x M" — >■ [0, oo) given by g{r, p) := |pp + ii'|rp. Since i is a diffeomorphism, 
this implies that Xh is also complete (see [19, Lemma 1.6.4]). 

(C) Many other examples with WH = g{H) can be obtained using homogeneous Hamiltonians functions. 
For instance, consider := q'^/q^ + q^/q'^ and ^{q,p) := piq^ +P2q^ on M := (M^ \ {0}) x M^. 

Then one has VH = H"^ - A, (pt{q,p) = {q,p- t^{q,p)) and 

Crit(ff) = Crit(F, $) = e \ {0} | q^ = ±q'^} x M^. 

4.2 H = h{p) 

Consider on M := R^" a purely kinetic Hamiltonian H{q,p) := h{p) with h € C°°(R"; R), and take the usual 
position functions ^{q,p) := q with d = n. Then ipt{q,p) = {q + t{'Vh){p),p), \/H = V/i, and Assumption 

3.2 is satisfied: 

{{<^,,H},H] ^ {{d,h){p)Mp)] 
In this example, we have Qr\\.{H) = Qnt{H, = R" x (V/i)-i({0}). 

4.3 The assumption {{^j,H},H} = as a differential equation 

Consider on M := R^" an Hamiltonian function H with partial derivatives Hp^ := dH/dpk and H^k := 
dUjdc^ . Then, finding the functions $j of Assumption 3.2 amounts to solving for $o the second-order linear 
equation 

As observed in Section 3.1, this is essentially equivalent (when k independent first integrals J\ = H, J2, ■ ■ ■ , Jk 
are known) to find the solutions $0 to 

n 

{$0, H} = Y, {Hp.d^e - Hgedp,)^o = 9{Ji, • . • , Jk)- (4.1) 

1=1 

The case g = 1 is sufficient, though trying to solve {$0; = 1 can at best provide solutions which are C°° 
outside the set Qr\t{H). A way to remove these singularities could be to multiply the solutions by a function 
g{H) that vanishes and is infinitely flat on Crit(iJ). For instance, if H[ Crit(i7)) consists of a finite number of 
values ci, . . . , Cg G R, one could take g{H) — n^=i e^(^^'^j) . Another possibility is to restrict the study to 
a submanifold M' of M (typically an open subset of the same dimension). However, problems can arise as the 
same (induced) symplectic structure (or Poisson bracket) must be used for the dynamic to remain unchanged; 
in particular, it must checked that the Hamiltonian flow preserves M' . 

(A) Repulsive harmonic potential. In this example we first solve the equation {$0, H} = 1, and then correct 
the functions $0 to make them C°°. So, let us consider fox K the Hamiltonian H{q,p) := ^ — 
K'^\q\'^) onM := R^". One can check that Crit(F) = {0} and that 

For j e {l,...,n}, take := ta,nh~^{Kqi /pj), where tanh"^(^) = ^Inji^j is C°° on 

R \ {±1}. Whenever pj — ±Kq\ the $j are not well-defined, but outside these regions, they satisfy 
{$j , H} = 1. It is possible in this case to get rid of the singular regions. Indeed, the functions Hj {q, p) := 
\ (p'j — K'^{q^Y) are first integrals of the motion and the singular regions correspond to the level sets 

Hj^{{(i}). Therefore, the functions := e~^3 i>j are well-defined and satisfy Assumption 3.2: 

{{$;,if},il} = {e-V,ij}=0. 
In this example, one has {0} = Crit(if) C Crit(i?, = fjj -^^^({O}). 
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(B) Simple pendulum. In this example we first consider the dynamics on a manifold and then restrict it to an 
appropriate submanifold. For K > 0, take H{q,p) := + K{1 — cosg)) on M := M?. One has 

Crit(-ff) = ttZ X {0} (the values q G I-kTL correspond to minima, while q E 27rZ + w correspond to 
inflexion points). Then, consider the open subset M' of M defined by the relation H > K, i.e. M' := 
{((7,p) e I — K cos^{q/2) > O}. One verifies easily that M' is preserved by the Hamiltonian 
flow, that M' n Crit(iJ) = and that M' corresponds to the region where the values of q along an orbit 
cover all of M. Define also 



where F{-\-) denotes the incomplete elliptic integral of the first kind. Then one verifies that the function 
$ is well-defined on M' and a direct calculation gives {$, H}{q,p) = p/\p\ for each {q,p) S M'. Now, 
p/\p\ = 1 on one connected component of M' and p/\p\ = — 1 on the other one. Thus Assumption 3.2 is 
verified on M' and Cr\t{H, $) = 0. 

(C) Unbounded trajectories of central force systems. Once again, we first consider the dynamics on a manifold 
and then restrict it to an appropriate submanifold. For K E M. \ {0}, take H{q,p) := ^(bP — ^l^l"^) 
on M := (M" \ {0}) x M", with n > 1 if if > and n > 1 if iiT < 0. One has Cnt{H) = 0. 

When K > (and n > 1), we must restrict our attention to the case where the Hamiltonian function 
H is positive (to avoid periodic orbits), and where at least one of the two-dimensional angular momenta 
Lij {q, p) := q'-pj — q^pi is nonzero (to avoid collisions, i.e. orbits whose flow is not defined for all f S M, 
see L29J). Therefore, the open set M' := {{q,p) G M\H{p,q) > 0, J2lj=i \Lij{q,p)\'^ ^ O} is an 
appropriate submanifold of M when K > 0. 

Consider now the real valued functions on M (resp. M') when K <0 (resp. K > and n > 1) given by 

^-^'^'^^ ■= dip) ^ 2{2Hiq,p)f^ {mHi,,P) + ± 2Vm^)P ■ . 
Since |pp < 2H{q,p) (resp. |pp > 2H{q,p)), then 

{^/2H{q,p)- \p\)^ > =^ 2H{q,p)\p\^ ± 2^2H(q,p)p ' > 



^ \q\{2H{q,p) + \p\^) ± 2^2H{q,p)p ■q>Q. 

So, $± are well-defined, and further calculations show that {$±, = 1 on M (resp. M'). As before, 
Crit(if) = Crit(if, $±) = 0. Note that $±(g,p) = p ■ g/bP when K = Q, which is coherent with the 
canonical function $ for the purely kinetic Hamiltonian H{q, p) = | 

One can construct a more intuitive function $o in terms of $±, namely, 



which also satisfies {<&o, H} = 1. Since the functions satisfying Assumption 3.2 are linear in t, one can 
regard them as inverse functions for the flow. The appearance of the inverse hyperbolic function tanh~^ 
in $0 is related to the fact that unbounded trajectories of the central force system given by if > are 
hyperbolas. 

(D) Poincare ball model. Consider i?i:={geK"||q'|<l} endowed with the Riemannian metric g given 
by 
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LetM := T*Bi ~ {{q,p) G Bi 
S^j=i d?-' A dpj, and let 



be the cotangent bundle on Bi with symplectic form 



H:M^R, {q,p) ^ \ ^ g^\q)PjPk = - 

be the kinetic energy Hamiltonian. It is known that the integral curves of the vector field Xh correspond 
to the geodesies curves of {Bi,g) (see [19, Thm. 1.6.3] or [11, Sec. 6.4]). Since, {Bi,g) is geodesicaUy 
complete (see Proposition 3.5 and Exercice 6.5 of [25]), this implies that Xjj is complete. There remains 
only to find a function $ satisfying Assumption 3.2 in order to apply the theory. 

Some calculations using spherical-type coordinates suggest the function 



$ : M 



Indeed, since 



(p-'z)(i-kP 



Vm<hp){i + \q\') 



2(P • q) 



bl(i 



< 



2|<7l 



l + \q\ 



<1, 



the function <1> is well-defined. Furthermore, direct calculations show that <& is C°° and that {<&, H} = 
e-^/" V2H. Therefore, Assumption 3.2 is verified and one has Crit(if) = Cnt{H, $) = Si x {0}. 

In one dimension, q{t) := tanh(f) is (up to speed and direction) the only geodesic curve, and 

^q,p) = e-V«(9.p) tanh-i ( , A = 2e-i/«(9.p) P tanh-i(g). 

Vbl(i + r)/ bl 



So, apart from the smoothing factor 2 e our $ coincides in one dimension with the inverse function 
of the flow. 



4.4 Passing to a covering manifold 

In this subsection we briefly discuss a way of avoiding the obstruction of periodic orbits: Given M a symplectic 
manifold with symplectic form co and Hamiltonian if, we let tt : M M \ Crit(il) be C°°-covering manifold. 
In order to preserve the dynamics, we endow the manifold M with the puUback u> := Tr*u) of the symplectic 
form uj and with the pullback i? := 7r*i? of the Hamiltonian H? 

Here are two simple examples of finite-dimensional symplectic covering manifolds. 

(A) Consider on the sphere M := (as seen in and with its standard symplectic structure) the Hamil- 
tonian H given by the projection onto the ^-coordinate. Outside the 2 polar critical points, all the orbits 
are periodic: the flow corresponds to rotations around the z-axis. In this case, one can use the covering of 
§2\{(0,0,±l)}givenbyM := {{'d,z) | i? G R, z G (-1, 1)} and the covering map 

TT : M M \ Crit(iJ) = \ {(o, 0, ±1)}, (i?, z) ^ (Vl - z'^ cos(i?), ^^l-z'^ sin(i?), z). 

Consequently, if : M ^ (— 1, 1) is the projection onto the z-coordinate and w = di? A dz. One can also 
check that ipt z) := (i? + 1, z) is the flow of H and that {$, H} = 1 for $(■!?, z) := So, Assumption 

3.2 is verified on M and Crit(_/7) = Cnt{H, $) = 0. 

'if one wants to consider only a Poisson manifold M, a Poisson structure can also be defined on M given that tt is C°^. Indeed, for 
U C M \ Crit(-ff) a sufficiently small open set (i.e. such that ■k~^(U) is a disjoint union of diffeomorphic copies), connected components 
of 7r~ ^ (J7) are diffeomorphic to U and the Poisson structure can be induced by this diffeomorphism. 
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(B) Harmonic oscillator. Consider on M := M^" (with its standard symplectic structure) the Hamiltonian 
given by if := ^{\p\'^+K'^\q\'^), where K e K\{0}. Define M := {(r,i?) | r e (0,oo)", i) e K"} 
and TT : M M \ Crit(iJ) = M^" \ {0}, with 

7r(r, t?) := {K~^ri cos(i?i), . . . , if" V„ cos('!?„), n sin(?9i), . . . , r„ sin(i?„)) . 

Then H{r,'d) = ijrp, w i^-^ X^jLi drj A d^j, and V5t(?-,i?) = (r,'d - Kt) is the flow of H. 
Furthermore, one has {<5?j, if} = —K for each function $j(r, t?) := dj. Therefore, Assumption 3.2 is 
verified on M with $ = ($i, ...,$„) and Crit(^) = Cnt{H, $) = 0. 

4.5 Infinite dimensional Hamiltonian systems 
4.5.1 Classical systems 

In the following examples, the infinite dimensional manifold M is either [? (K) or (R) © (K) (equivalence 
classes of real valued square integrable functions)'*. The atlas of M consists in only one chart, the tangent space 
T„M at a point u e M is isomorphic to M, and the Riemannian metric on M is flat {i.e. independent of the 
base point in M) and given by the usual scalar product ( • , • ) on (R) or \} (R) ® (K) . 

To define the symplectic form on M in terms of the metric ( • , • ) we let "H", s S R, denote the real Sobolev 
space or ® (see [4, Sec. 4.1] for the definition in the complex case) and we let y denote 

the real Schwartz space ,S^{K) or ,i^(R) © ,S^{K). Then we consider an operator J : S/ .5^ (which can be 
interpreted by continuity as an endomorphism of the tangent spaces TuM ~ M) satisfying the following: 

(i) There exists a number dj > 0, called the order of J, such that for each s G R the operator J extends to 
an isomorphism 'H^~'^-' (which we denote by the same symbol). 

(ii) J is antisymmetric on i.e. {Jf, g) = — (/, Jg) for all f,g€y. 

It is known [22, Lemma l.lj that the operator J := —J~^ : M 'H'^^ (of order —dj) is bounded and 
anti- self adjoint in M. In consequence, for each s > the map w : 'H* x 'H* ^ R given by 

w(/,<?) :=-(J/,<?) 

defines a symplectic form on 'H*. 

The functions on the phase space (such as H or $j) are infinitely Frechet differentiable mappings from 

(a subset of T-L^" for some sh > 0) to R, i.e. elements of C°°{Osh ! The Hamiltonian function H S 
C°°(e's^;R) is defined as follows: for some h G C°°(R*^+i;R) (or h G C°°(R2('=+i);R) if M = L^{R) © 
L2(R)), one has for each u € Osh 

H{u):= I dxh{uQ,ui,...,UsH), 
Jr 

where Uj := ^ . Since H G C°°(e's„;R), the differential of if at m G Og^ on a tangent vector f e y c 
M ~ T„M is given by 

t^ot^ J Jr ^ duj dx^ ^ Jr dx^ duj 

where the second equality is obtained using integrations by parts (with vanishing boundary contributions). 
The (Riemannian) gradient vector field gradil associated to the linear functional dH satisfies by definition 
((gradi?)(u), /) = dHu{f) for all u e Osh f €: y (here {gradH){u) a priori only belongs to the 

''in the case of the wave and the Schrodinger equations below, one can easily extend the results to the situation where L^(M) is replaced 
by L'^(IR"). We restrict ourselves to the case n = 1 for the sake of notational simplicity. 
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topological dual y* of 5^ , which means that ( • , • ) denotes a priori the duality map between y* and S^). So, 
(grad-ff)('u) is given by 

(grad/.)(.)= £(-1)^1-1^. (4.2) 

Then, the Hamiltonian vector field Xn is the map ^ -5^* satisfying 

(J/,Xh(«)) = = dHu{f) = (/, (gradif)(«)) 

for all u G and / e J^. Since J is anti- self adjoint, this implies that JXh{u) = — (gradil)(w) in 
which is equivalent to Xh{u) = J{gradH){u) in S^* . So, the equation of motion with Hamiltonian 

H has the form ^ u = J{gra,dH){u), and {^,H} = di>(X^f) = (grad$, J(gradi?)) for all functions 
e C°°{Osu ; K) with appropriate gradient. 

Before passing to concrete examples, we refer to [20] for standard results on the local existence in time of 
Hamiltonian flows (global existence is specific to the system considered). 

(A) The wave equation. We refer to [1, Ex. 5.5.11, [2, Ex. 8.1.12], [13, Sec. 2.1] and [30, Sec. X.13] for a 
description of the model. The existence of the flow for all times depends on the nonlinear term in the 
Hamiltonian (see for instance [30, Thm. X.74] and the corollary that follows). 

In this example, the scale {'W'}s>o is given by := ?^*(M) © n^R). The metric on M := L^(R) ® 
L^(IR) is given for each (p, q), {p,q) €. M by ((p, q), {p, q)) := da; {pp + qq), and the operator J is 
given by 

J:M^M, ip,q)^i-q,p). 

It is an isomorphism of degree with J = J. Given m > and F g C°°(]R; M), one can find a subset 
Oi CH^ (depending on F) such that the Hamiltonian function 

H-.Oi^R, {p, q)^ [ dx h{p, q, d^q) = \ [ dx {p^ + [d^qf + m^q^ + 2F(g)}, 

JR JR 

is well-defined and C°°. In fact, we assume that 0\ is chosen such that (i) all the functions on the phase 
space appearing below are elements of C°° (d ; K), and (ii) integrations by parts involving these functions 
come vanishing boundary contributions. Then one checks that (gradi? ) (p, q) ~ (p, rr?q + -F'(g) — 
due to (4.2), and that Xh{p, q) is trivial if and only ifp = Q and m?q + F'{q) — d^q = 0. The constraint 
on q depends on the choice of F. For example, when F{q) = 0, g or q^, the solution q of the differential 
equation does not decay as ^ oc . In consequence, the corresponding pairs (p, q) cannot belong to M, 
and Crit(if) = {(0, 0)}. The equation of motion 

^(P,<?) = J(gradF)(p,9) (4.3) 

coincides with the usual the wave equation since the combination of ^p = d^q — m^q — F'{q) and 
^q=p gives 

q = d^q - m^q - F' (q). 

When m ^ 0, this equation is called the Klein-Gordon equation, and F is usually assumed to be a 
nonUnear term of the form F{q) = q^ for some A e M. A first relevant observation is that the function 

Co G C°°(C'i; M) given by Co(p, q) := dxp{dxq) is a first integral of the motion. Furthermore, the 
function $0 e C°°(Oi; M) given by ^o{p, q) := dx id^ h{p, q, d^q) has gradient (grad<5?o)(p, l) = 
( idR p, idg rri^q + id«. F'{q) - dx{idm d^q)) . Therefore, 

{^Q,H]{p,q) = ((grad$o)(p,9),^(gradi?)(p,g)) = / dx p [id^ dlq - dx{\dB.d^q)] = -Co{p,q), 
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and $0 satisfies Assumption 3.2. Here, we clearly have 

Crit(i?,$o) = Co\{0}) = {{p,q) e Oi I J^p{d.q)dx = O} D {(0,0)} = Crit(ff). 

If we assume further that F = 0, then the equation of motion (4.3) is linear. Therefore any pair (9|p, dlq), 
j > 1, with {p, q) a solution of (4.3), also satisfies (4.3). Consequently, if the subsets Oj C W have 
properties similar to the ones of Oi, then the functions Cj € C°°(C'j+i;R) and Hj € C°°(Oj+i;M) 
given by Cj (p, g) /jj da; (S^p) (9^+^g) and Hj{p,q) := J^dx h(^dlp,diq,di'^^q) are first integrals 
of the motion. Accordingly, one deduces that the functions <J>j € C°°{Oj+i;R) given by ^j{p, q) := 
Jjjda; idm. h[dlp,diq,di'^^q) satisfy {^j,H} = —Cj on Oj+i. So, if F = 0, there is an infinite 
family of functions satisfying Assumption 3.2, and one has again Cr\t{H, $j) D Crit(i?), with 
dl : Crit(iJ, $j) ^ Crit(iJ, $o) an isomorphism. 

Finally, when F = and m = one can check that the function $0 G C°^(C'i; M) given by ^o{p, q) '■= 
da; ids, p{dxq) has gradient (grad$o) {p, q) = (idiR d^q, - ids O^p - p). Then, 

{$0,i?}(p,g)= / dx{id^{d^q){dlq)-idRpd^p-p')=-\ I dx{{d,qf+p^)=-H{p,q), 
Jr Jm. 

where the third equality is obtained using integrations by parts (with vanishing boundary contributions). 
Thus $0 satisfies Assumption 3.2. Furthermore, since {^0, H}{p, q) = implies dx { {dxqY' = 
0, one has Crit(iJ, $0) = Crit(i?) = {(0,0)}. As before, any derivative of a solution of the equa- 
tion of motion is still a solution of the equation of motion. So, it can be checked that the functions 

e C°°(0,+i;K) given by ^j{p,q) ^dx idg {dip){di+\) satisfy {^j,H} = -Hj on Oj+i. 
Therefore, one has once again Crit(i7, $j) = Crit(7J) = {(0,0)} and the $j's constitutes a second 
infinite family of functions satisfying Assumption 3.2. 

(B) The nonlinear Schrodinger equation. We refer to [22, Ex. 1.3, p. 3 &5] for a description of the model. 
The existence of the flow for all times depends on the nonlinear term in the Hamiltonian (see for instance 
[8, Sec. 1.2] and [33, Sec. 3.2.2-3.2.3]). 

The setting is the same as that of the previous example, except that the Hamiltonian function H e 

C°°(e'i;M) is given by 

H{p, q):=\[ dx{{d^pf + {d^qf + V-{p'' + q^) + F{p^ + q^)}, 

where V,FeC°° (R; R). Using (4.2), one checks that the gradient of H at {p, q) e Oi is 

{gradH){p, q) = {-dlp + Vp+pF'{p^ + q^), -d^q + Vq + qF'{p^ + q^)). 
So, the equation of motion ^{p,q) = J(gradi? ) {p, q) is equivalent to the nonlinear Schrodinger equation 

^u = i{-dlu + Vu + uF'{\u\^)), (4.4) 

with u := p + iq. Without additional assumptions on F or V, it is hardly possible to determine the set 
Crit(iJ) of functions u for which the r.h.s. of (4.4) vanishes. However, it is known that in general Cnt{H) 
is not trivial, as in the case of eUiptic stationary nonlinear Schrodinger equations (see Theorem 1.1 and 
Proposition 1.1 of [6]). 

Now, assume that V = F = and for each j > 1 let Oj C W be a subset having properties similar 
to the ones of d. Then the functions Hj e C°°{Oj;R) and Cj G C°^{Oj+i;M.) given by Hj{p, q) := 
U^dx{idiqr + idip)^} ^ J^dxh,ip,q)mdC,ip,q) := J^dx {{d^qm+'p) ~ idt'q)idipl} = 
/ju da; Cj {p, q) are first integrals of the motion. Furthermore, the functions £ C°° {Oj ; R) and <&j G 
C°° {Oj+i ; M) given by {p, q) := dx ida hj (p, q) and $j (p, q) := da; idg Cj (p, q) satisfy {^j,H} = 
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Cj and {^j,H} = 4ffj+i on O^+i. So, the $j's and the $j's constitute two infinite families of func- 
tions satisfying Assumption 3.2. Note that the sets Cr\X.{H,^j) = Cj^{{lS\) = {{p,q) e Oj+i | 
Jgdx = 0} (with isomorphisms 9^ : Crit(if, $j) Crit(/f, $0)) are rather large, 

whereas Crit(i7, $j) = Crit(i7) = {(0,0)}. 

Some of the above functions still work when V and F are not trivial. For instance, the identity {<i>o, H} = 
Co on 0\ remains vaUd for all V and F. Furthermore, if y = Const., then {Co, H} = on Oi. Conse- 
quently, <l>o satisfies Assumption 3.2 for all F and for V = Const., and one has Crit(i7, $0) 2 Crit(i?). 
This last example is interesting since it applies to a large class of nonlinear Schrodinger equations. 

(C) The Korteweg-de Vries equation. Among many other possible references, we mention [1, Ex. 5.5.7] and 
[22, Ex. 1 .4, p. 3 & 5]. For the global existence of the flow, we refer the reader to [12, Sec. 1] and references 
therein. 

In this example, the scale {'H^}s>o is given by := 'W(M) and the sets Oj, j G N, are appropriate 
subsets of W . The Hamiltonian function € C°° (d ; K) is given by 

H{u) := f dx{^{d^uf + u''), 

and the isomorphism J := 9^ is of order 1. 

The gradient of at u G Oi is —d'^u + 3w^. So, the elements of Crit(_ff) are functions u satis- 
fying —d^u + 3u^ = 0; these are Weierstrass p-functions [21, Sec. 134. F], that is, functions with 
many singularities and no decay at infinity. Thus, Crit(fl') = {0}. Furthermore, the equation of motion 
= J{gra.dH){u) coincides with the KdV equation jj-u = dx{ — d^u + Su^) . 

There exists an infinite number of first integrals of the motion with polynomial density, that is, of the form 
Hj da; hj, where hj is a finite polynomial in u and its derivatives (see [27, Sec. 3]). For example, 

when hi{u) = u, h2{u) = it^, hz{ti) = ^{d^uf + u^, or h4{u) = {d'^uf + 10u{dxuf + Su". So, let 
$0 G C°°(Oo;R) begivenby $o(w) := Js_'ix idg u. Then the gradient of $0 at wis idR, and {$0, -f^} = 
— 3i?2 on Oi. Since H2 is a first integral of the motion, this implies that $0 satisfies Assumption 3.2. 
Furthermore, the fact that H2{u) = ||u||l2(m) implies that Crit(il, $0) = {0} = Cnt{H). 

Looking for others $ of the form = da; g(x) G(u, dxU, . . . , d^u), with G a polynomial and g 
a C°° function, is unnecessay. Indeed, both {$, H} and T{t) := $ — t{$, H} are first integrals of the 
motion with density C°° in x and polynomial in u and its derivatives (and t-linear in the case of T). Thus, 
we know from [34, Thm. 1 & Rem. 3] that they are completely characterised, up to the usual equivalence 
of conservation laws [28, Sec. 4.3]. Therefore, the functions $ are also completely characterised. Note 
however, that it is not excluded that functions $ with an integrand G involving fractional derivatives, 
an infinite number of derivatives, or of class might work. Non-polynomial conserved densities are 
known to exist in the periodic case (see [27, Sec. 5]). 

4.5.2 Quantum systems 

Let be a complex Hilbert space, with scalar product ( • , • ) antiUnear in the left entry. Define on H the usual 
quantum-mechanical symplectic form 

uj-.HxH^R, (V'i,V2) 2lm(V'i,V2). 

The pair {H,oj) has the structure of an (infinite-dimensional) symplectic vector space. Now, define for any 
bounded selfadjoint operator Hop G ^{H) the expectation value Hamiltonian function 

Then, it is known [26, Cor. 2.5.2] that the vector field and the flow associated to H are Xh = —iH^p and 
</'t(V') = e~'*^°p tl). Therefore, the Poisson bracket of two such Hamiltonian functions H, K satisfies for each 

{K,H}{'ip) = io{XK{il'),XH{''P)) = -uj{Kop^,Hopi>) = {^,i[Kop,Hop]^). 
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So, in this framework, verifying Assumption 3.2 amounts to finding Hamiltonian functions H = {Hop) and 
$j = (($j)op) satisfying the commutation relation 

)op, ifop], ^?op] =0. (4.5) 

In concrete examples, the operators Hop and (<&j )op are usually unbounded. Therefore, the preceding calcula- 
tions can only be justified (using the theory of sesquilinear forms) on subspaces of H where all the operators 
are well-defined. We do not present here the whole theory since much of it, examples included, is similar to that 
of [31]. We prefer to present a new example inspired by [16], where all the calculations can be easily justified. 

Let U be an isometry in H admitting a number operator, that is, a selfadjoint operator A'' such that UNU* = 
N — 1. Define on H the bounded selfadjoint operators 

A := Re(C/) = and S := \m{U) = ^^{U - U*). 

Then we know from [16, Sec. 3.1] that any polynomial in U and U* leaves inviiriant the domain ViN) C Hof 
N. In particular, the operator 

Ao:=^{SN + NS), V{Ao) := V{N), 

is well-defined and symmetric. In fact, it is shown that Aq admits a selfadjoint extension A with domain 2?(^) = 
2?(-/VS'). Furthermore, one has on 'D{N) the identity «[j4, A] = A^ — 1. So, if we define the Hamiltonian 
functions 

H-.n^R, V^(A)(V') and <i>:V{N)^R, ip ^ {A){il)), 
we obtain for each tp G 2?(iV) 

{Wnm = = l<i[A,H]){i,) = (A2 - 

and Assumption 3.2 is verified for each ij) e 'D{N): 

{{$, H}, H}{^) = a;(X<A2_i)(V), = (^[A^ - 1, A])(V) = 0. 

Now, since the spectrum of A is [—1,1], the operator 1 — A^ is positive, so we have the equivalences 

(A2-1)(V')=0 <^ 11(1- A2)i/2i/,f =0 ^ ij€E^{{±l}). 

Thus, 

Cnt{H, $) = {WH)-\{0}) = e V{N) \ (A^ ^ = O} = ©(A^) n E^{{±1}). 

On the other hand, the elements tjj £ Crit(if) satisfy the condition 

= Xh{iP) = -iAip <^ Ve£;^({0}). 

This implies that Cnt{H) = {0}, since the spectrum of A is purely absolutely continuous outside the points ±1 
[16, Prop. 3.2]. Finally, the function Tf is given by 

T/ = -(^).(Vi?/)((A2-l)) 

onr>(A^) \ Crit(iJ,$). 

Typical examples of operators A and N of the preceding type are Laplacians and number operators on 
trees or complete Fock spaces (see [16] for details). 
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